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^ ! Accurate astrometry and photometry of saturated and coronagraphic point 

\0 '. spread functions (PSFs) are fundamental to both ground- and space-based high 



contrast imaging projects. For ground-based adaptive optics imaging, differential 
atmospheric refraction and flexure introduce a small drift of the PSF with time, 
and seeing and sky transmission variations modify the PSF flux distribution. 
For space-based imaging, vibrations, thermal fluctuations and pointing jitters 
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c^ ! can modify the PSF core position and flux. These effects need to be corrected to 

properly combine the images and obtain optimal signal-to-noise ratios, accurate 
k> i relative astrometry and photometry of detected objects as well as precise detec- 

H ! tion limits. Usually, one can easily correct for these effects by using the PSF core, 

but this is impossible when high dynamic range observing techniques are used, 
like coronagrahy with a non-transmissive occulting mask, or if the stellar PSF 
core is saturated. We present a new technique that can solve these issues by using 
off-axis sateUite PSFs produced by a periodic amplitude or phase mask conju- 
gated to a pupil plane. It will be shown that these satellite PSFs track precisely 
the PSF position, its Strehl ratio and its intensity and can thus be used to register 
and to flux normalize the PSF. A laboratory experiment is also presented to val- 
idate the theory. This approach can be easily implemented in existing adaptive 
optics instruments and should be considered for future extreme adaptive optics 
coronagraph instruments and in high-contrast imaging space observatories. 

Subject headings: Instrumentation: AO - planetary systems - stars: imaging 
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1. Introduction 

Accurate astrometry and photometry of saturated and coronagraphic point spread func- 
tions (PSFs) are fundamental to botli ground- and space-based liigli contrast imaging of 
brown dwarfs and exoplanets. Achieving precise astrometry is important to minimize the 
time required to confirm via proper motion analysis that two objects are gravitationally 
bound, i.e. that a newly found companion is not a background object, and to better con- 
strain its orbit. Precise photometry is essential to better understand the companion physical 
characteristics. For ground-based adaptive optics imaging, both standard observing tech- 
niques and more specialized observing methods, like simultaneous spectral differential imag- 
ing (Racine et al. 1999; Marois et al. 2000; Sparks & Ford 2002; Biller et al. 2004; Marois 
2004; Marois et al. 2005), angular differential imaging (Marois 2004; Liu 2004; Marois et al. 
2006) and coronagraphy (Lyot 1932), generally require the acquisition of a long sequence of 
images. The object airmass, the telescope orientation and the observing conditions (seeing 
and sky transmission) change during the sequence, producing point spread function (PSF) 
flux and Strehl ratio variations over time. The PSF is also slowly moving and evolving 
due to flexure and differential atmospheric refraction if a near-infrared camera is used while 
operating the adaptive optics wavefront sensor in the visible. Space-based imaging is also 
affected by similar problems from vibrations, thermal fluctuations and pointing jitters. If 
not corrected, these effects bias the relative photometry and astrometry of detected sources 
and the estimated sensitivity limits. To obtain the optimal signal-to-noise ratio (S/N) and 
precise astrometry and photometry after combination of all acquired images, each individual 
image needs to be accurately registered and flux normalized. 

Standard techniques to properly register a PSF and to estimate the relative astrometry 
and photometry generally require the PSF core. If the PSF core is saturated to improve 
observing efficiency or if a non-transmissive focal plane occulter is used, the PSF core is 
not visible and other techniques need to be applied. If available, an off-axis sharp ghost 
image can be used to register the PSF and estimate its core position, but the possibility of 
a non-common motion between the ghost and the PSF makes the registering uncertain. 

The usual technique for estimating sensitivity limits (or a source relative intensity) when 
the PSF core is saturated or occulted is to compare the noise in the combined image (or 
a source peak intensity) to the peak intensity of an unsaturated/unocculted PSF acquired 
before and/or after the saturated/occulted image sequence. Such estimates can be biased 
by the short exposure time and the generally small number of unsaturated/unocculted im- 
ages. Furthermore, the observing conditions of the unsaturated/unocculted images can be 
signiflcantly different from those of the saturated/occulted image sequence. 

To obtain the highest S/N ratio in the combined image, one would want, ideally, to 
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normalize the images such that any point source, after image registration, has the same 
number of counts per full width at half maximum (FWHM) in each image. The images 
would then be weighted by the inverse square of their noise and combined. 

A simple solution to these astrometric and photometric problems exists. A periodic 
phase or amplitude mask can be introduced at a pupil plane to produce fainter off-axis 
copies of the primary PSF. The position and relative intensity of these satellite PSFs are 
mainly fixed by the mask amplitude and spatial periodicity. These satellite PSFs track 
precisely the PSF core position as well as its intensity variations. They can thus be used for 
accurate astrometry and photometry. First, this technique will be analyzed using a simple 
analytical model and numerical simulations. The technique will then be validated using a 
laboratory experiment. A technique similar to the one presented in this paper has been 
independently developed by Sivaramakrishnan & Oppenheimer (2006). 



2. Non-Coronagraphic PSFs 

The effect of a pupil-plane conjugated aberration can be estimated using a far- field 
approximation. The intensity /(r^,^), where r] and ^ are coordinates in the focal plane, of 
the PSF is simply the Fourier transform FT of the complex electric field (Schroeder 1987) 
at the pupil 

/(r^,0 = |FT(A(x,y)e^[^(^'^)l)f = \F,f ^ h (1) 

where A is the wavefront amplitude, equal to one inside the pupil and zero elsewhere, and 
its phase, having a noise RMS a inside the pupil. Both A and are expressed in coordinates 
X and y in the pupil and are real functions. The phase has zero mean. The functions 
la and Fa are respectively the aberrated PSF intensity and its complex focal plane electric 
field. The next two sub-sections will study the effect of adding a specific phase or amplitude 
pupil-plane mask to the aberrated wavefront. It is assumed in the following analysis that 
(t>{x,y) is small (space-based or AO imaging), producing diffraction limited images. 



2.1. Pupil-Plane Periodic Phase Mask 

If a phase aberration 0mask having a FT equal to $mask and a standard deviation o"mask 
is introduced on the wavefront at a pupil plane, the modified PSF intensity can be deduced 
from Eq. 1. The aberration is simply replaced by + 0mask and the following expression 
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I{r], = |FT(A(x, y)ei[^(^'y)+<^--k(x,y)]) |2 _ ^2) 

Separating the effect of the mask, we find 

7(77,0 = |FT([A(x,y)e'['^("'y)l][e'['^— ^("'y)l])f . (3) 

Following the work of Bloemhof et al. (2001); Sivaramakrishnan et al. (2002); Perrin et al. 
(2003), assuming 0mask is small, the complex exponential of 0mask can be expanded using a 
Taylor approximation to find 
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FT([A(x,y)ei[^(^'y)l][l + i0„,,3k - ^ + ..•]) 



We can then use the equation derived in Perrin et al. (2003) to find 
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Hv,0 = J2p- (5) 



n=0 
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P« = ^" E kUnLkvS ^^ ^' <^mask)(F: ^"-'^ $Lsk)- (6) 

fc=0 '^ ^' 

If the phase mask is chosen to be a periodic function, $mask is sharply peaked at symmetric 
locations. The terms proportional to Fa ** $mask thus simply add replicas of Fa at each peak 
of $mask and at integer multiples of these peaks for multiple convolutions. Keeping only 
terms up to the second order, we find 



1(7], 0=4 + 23[F;(F, ^ $^ask)] - 5R[F;(F, ^ $,nask ^ $mask)] + \Fa ^ ^maskl' (7) 

where the symbols Q" and K denote respectively the imaginary and real parts. Both the 
second and third terms are interference effects between the aberrated PSF and $mask- For 
the replicas to be greater than la at the same separation, o"mask needs to be larger than 
the noise component of a at the spatial frequency corresponding to the mask periodicity. 
Nothing that Fa sharply drops with angular separation, since both the second and third 
terms are multiplied by Fa and since $mask is selected to be brighter than Fa ($mask ^ -^a 



-5- 



at $mask maxima), both terms can be neglected when compared to the fourth one. In that 
case, we find that the PSF intensity evolution due to the phase mask is simply 

/(^,0-4=l^a^$mask|'. (8) 

Since Fa is the aberrated PSF electric field, the periodic phase mask produces replicas of 
the aberrated PSF, we call these satellite PSFs. Since the satellite PSFs are produced by a 
convolution of the aberrated PSF, they track the PSF core position and intensity variations. 
The intensity and location of the satellite PSFs are respectively set only by the amplitude 
and periodicity of the phase mask, they are thus insensitive, to first order, to a mask offset. 
A discussion on the neglected terms and their effects on the astrometric and photometric 
precision can be found in sect. 4. 

If we assume that the phase mask is a sine wave, two satellite PSFs will be produced 
at an angular separation equal to the number of cycles per pupil in X/D units and with 
a relative intensity equal to fT^iask/2- Fig. 1 shows an unaberrated 1.6 /xm monochromatic 
PSF with a 11 nm amplitude sine wave phase mask having 25 cycles per pupil as well as a 
PSF having the same sine wave phase mask but with a 150 nm RMS phase aberration. The 
aberration power spectral density (PSD) follows a power-law of index -2.6. It is expected 
that satellite PSFs will be 1000 times fainter than the PSF and located at 25 X/D. PSFs are 
simulated with a 256 pixel diameter circular pupil inside a 1024 x 1024 pixel image producing 
a PSF having 4 pixels per X/D. The satellite PSFs are clearly at their predicted location, 
they have the estimated relative intensity and they track the PSF core position, Strehl ratio 
and its intensity. These satellite PSFs can thus be used for registration and photometric 
calibration (see sect. 5). 



2.2. Pupil-Plane Periodic Amplitude Mask 

A similar analysis can be done for a pupil plane conjugated amplitude mask. Introducing 
an amplitude mask emask having a FT equal to -Emask and a standard deviation o"mask on the 
wavefront conjugated to the pupil plane modifies Eq. 1 as follows 

l{r],0 = |TF([A + Ae]e^[^(^'y)l)p. (9) 

Reorganizing this equation, we find 
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Knowing that the FT of a multiphcation is the convolution of the two FTs, we find 

I{r],0 = Ia + 2?ii[Fa{F:^E*)] + \Fa^E\\ (11) 

Following the argument of section 2.1, the second term can be neglected and we find the 
modified PSF intensity to be 

I{r],0-Ia=\Fa*E\^ (12) 

The satellite PSFs are thus again copies of the primary aberrated PSF. See sect. 4 for 
a discussion about the neglected term and its effect on the astrometric and photometric 
precision. 

If the amplitude mask is a sine wave, the satellite PSFs have a relative intensity equal 
to o"^ask/2 and their separation is equal to the number of cycles per pupil in \/D units. 
Fig. 2 shows an unaberrated 1.6 ^m monochromatic PSF with a 4% amplitude sine wave 
transmission mask having 25 cycles per pupil as well as an aberrated PSF with the same 
sine wave mask but with 150 nm RMS phase aberration as in section 2.1. It is expected that 
satellite PSFs will be at 25 X/D separation and 1000 times fainter than the PSF. Similarly to 
what has been found for a phase mask, the satellite PSFs are at their predicted location, they 
have the estimated relative intensity and they track the PSF core position and its intensity. 

Instead of a sine wave amplitude mask, another option is to simply use a transmission 
amplitude grating consisting of wires to mask a periodic section of the pupil. In that case, 
Eq. 12 is still valid and the satellite PSFs will be the convolution of the wire grid Fourier 
transform with the aberrated PSF. The wire grid can be treated as a multi-slit Young 
experiment and the satellite PSFs intensity function /* is thus simply (Hecht 1998) 

m = is('^y('^' (13) 

\ p J \ sma / 

where 9 is the angular separation along the satellite PSFs axis, Jq is the intensity at ^ = 0, A^ 
is the number of wires across the pupil and (3 and a are respectively equal to {kb/2) sm6 and 
{ka/2) sin 9. The constant k is the propagation number (27r/A) and the values a and b are 
respectively the wire spacing and thickness. The first parenthesis in Eq. 13 is the intensity 
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envelope produced by the diffraction of a single wire (sine function liaving ~ A/6 FWHM) 
wliile tlie second parentliesis is tlie intensity modulation from the interference of diffracted 
light from all wires. Satellite PSFs (local maxima) appear at mX/a separations, where m is 
an integer. In contrast to the phase and amplitude sine wave masks that show two satellite 
PSFs, the transmission amplitude grating produces numerous satellite PSFs. The relative 
intensity, at ^ = 0, of the satellite PSFs is equal to the square of the ratio of the area masked 
by the grating over that in the unobscured pupil region 

^0 \Jp^^,,{A + Ae)dxdyJ 

where Jq is the aberrated PSF peak intensity. For a uniform pupil having numerous wires, 
the circular shape of the pupil can be neglected and the intensity ratio is approximated by 

If a ^ b, the ratio is simply {h/af' . 

Fig. 3 shows an unaberrated 1.6 /im monochromatic PSF simulation with a grating 
having 25 wires across the pupil, each of width equal to 1/250 of the pupil diameter, as 
well as an aberrated PSF (150 nm RMS phase aberration) with the same grating (for these 
simulations, the pupil diameter is 250 pixels and the image has 1000x1000 pixels). It is 
expected that the first satellite PSFs will be at 25 \/ D separation and 81 times fainter than 
the PSF. The added satellite PSFs are at their predicted location, they have the estimated 
relative intensity and they again track the PSF position, Strehl ratio and its intensity (see 
Sect. 5). 



3. Coronagraphic PSFs with added phase/amplitude mask 

The analysis presented in Sect. 2 was done using a far field approximation without the 
use of a coronagraph. If a coronagraph is used, the amplitude/phase mask needs to be 
introduced before the focal plane occulter so that the satellite PSFs be replicas of the unoc- 
culted PSF. If we consider a simple pupil plane Gaussian apodizer, it is easy to understand 
from Eq. 8 and Eq. 12 that the conclusions of Sect. 2 are still valid for a coronagraph since 
the PSFs are simply obtained with A being the apodized pupil instead of a uniform pupil. 
The apodized Gaussian coronagraph was chosen here since the PSF diffraction core is still 



detectable and can thus be used to verify that satelhte PSFs track the PSF core position 
and its intensity. The resuh would be the same for a band-limited (Kuchner & Traub 2002) 
or a standard Lyot coronagraph since the grating is introduced in a pupil plane before the 
focal plane mask occulter and because, similarly to the apodized pupil coronagraph, these 
coronagraphs produce PSFs which are dominated by the second order halo term (Perrin et 
al. 2003; Bloemhof 2004). A Gaussian apodized pupil for a 1.6 /xm monochromatic PSF was 
simulated using the transmissive amplitude grating with and without 150 nm RMS phase 
aberration (see Fig. 4). The pupil apodization was simulated by convolving the uniform pupil 
by a Gaussian having a FWHM equal to one quarter of the pupil diameter. The simulation 
confirms that the satellite PSFs do again track the PSF core position and intensity. 



4. Expected Astrometric and Photometric precision 

The aberrated PSF, the neglected terms of Eq. 7 and 11 and the rest of the terms 
in Eq. 6 show symmetries or antisymmetries that can, in the presence of aberrations or 
mask shifts, slightly affect the PSF center and intensity derived from the satellite PSFs. 
The task of finding the dominant terms in different regimes is complicated by the fact that 
Fa is complex; the imaginary and real parts of terms in Eq. 7 and 11 are not necessarily 
respectively antisymmetric and symmetric as in the unaberrated case. 

For simplicity let us first consider the case without aberrations. In this case. Fa is the 
field strength of a perfect PSF; it is thus real and symmetric. The first neglected term of Eq. 7 
is antisymmetric and affects the astrometric precision and the second one is symmetric and 
affects the photometric precision (higher order neglected terms from the Taylor expansion 
have a negligible effects since they are at least $mask times smaller than these two terms). The 
amplitudes of these terms are function of $mask- If $mask is antisymmetric, the second term 
is zero and the first term is maximal (astrometric error only), while if $mask is symmetric, the 
first term is zero and the second term is maximal (photometric error only). The symmetry 
of the mask is determined by its position in the pupil and may vary due to fiexure or 
removal/reinsertion. The astrometric error is proportional to the ratio of the amplitude of 
the first neglected term to the fourth term of Eq. 7, or -Fa/'^mask- The photometric error is 
proportional to the ratio of the amplitude of the second neglected term to the fourth term, 
or Fa- Thus astrometric measurements with brighter satellite PSFs are less affected by a 
mask offset. 

For the amplitude mask without aberration, the neglected term of Eq. 11 produces only 
a symmetric structure (photometric error only). If the mask is antisymmetric, the term is 
zero, the satellite PSF intensity is thus equal to the expected value, while if the mask is 
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symmetric, the term is maximal and a photometric error is present. Again, a mask offset 
can induce symmetry variations and cause a photometric error. The photometric precision is 
proportional to the ratio of the neglected term to the third term of Eq. 11, or Fa/ ^mask- The 
photometric precision is less sensitive to a mask offset when the satellite PSFs are brighter. 

As mentioned previously, since Fa is complex, the case with aberrations is much more 
complicated. Astrometric and photometric errors can arise from both imaginary and real 
parts. Furthermore, the residual symmetric and antisymmetric speckle noise from the PSF 
will also influence the level of accuracy. Intuitively, it is expected that brighter satellite PSFs 
will result in a better astrometric and photometric precision. 



5. Simulated Performances 

In previous sections, it was shown that periodic amplitude or phase masks can be used 
for astrometry and photometry calibration of saturated and coronagraphic PSFs. In realistic 
conditions, the satellite PSFs are overlaid over a background of speckle noise that can affect 
the determination of the satellite PSFs centers and intensities. It is thus important to adjust 
the intensity of the satellites PSFs to ensure that the desired accuracy is realized. Since we 
have assumed that the speckle noise is limiting the astrometric and photometric precision, the 
sky background, detector read and photon noises are neglected. Numerical simulations are 
now used to analyze the registering and photometric precision of the technique. Each satellite 
PSF is cross-correlated with a Gaussian to determine its precise location. The average X 
and Y positions of two symmetric satellite PSFs are then used to derive the PSF center. 
This derived position is then compared to the calculated PSF center using the same cross- 
correlation algorithm on the PSF core. Relative intensities are simply obtained by taking 
the ratio of the satellite's peak intensity over that of the primary PSF. The coronagraph 
was simulated by simply convolving the pupil by a Gaussian having a FWHM equal to one 
quarter of the pupil diameter. The brightness of the satellite PSFs was adjusted to 0.1% of 
the PSF peak intensity, which correspond to 100 times the background speckle noise of an 
aberrated (150 nm RMS) non-coronagraphic PSF. Both the phase and amplitude sine phase 
masks are antisymmetric on the pupil while the grating is symmetric (it is impossible to get 
a purely antisymmetric mask using a transmission amplitude grating). The first simulation 
is performed without aberration (see Table 1). As mentioned in sect. 4, only the phase 
mask shows astrometric errors when the mask is displaced while all mask types show small 
photometric errors. The use of a coronagraph attenuates interference terms and yields a 
better astrometric and/or photometric precision in all cases. 

Twenty simulations using independent 150 nm RMS phase aberrations were then real- 



-10- 



ized for phase sine, amplitude sine and amplitude grating masks for both cases with and 
without the coronagraph. Table 2 summarizes the achievable registering and photometric 
precision. Typically, an astrometric error of 0.3 pixels RMS (~1/15 \/ D) is found for both 
the non- coronagraph and the coronagraph cases. Typical relative intensities are precise to 
less than 5%, enough for accurate flux normalization. 

Table 3 presents the registering and photometric accuracies for the non-coronagraph 
and coronagraph cases using the sine wave amplitude mask for different satellite PSF inten- 
sities relative to the local speckle noise. Satellite PSFs ~10 times brighter than the speckle 
noise limit the PSF registering accuracy to 1.2 pixels RMS (1/3 \/D) and the photometric 
calibration to 15-20%, while satellite PSFs being ~100 times brigther than the background 
speckle noise limit the PSF registering to 0.3 pixel RMS (~1/15 A/D) and the photometric 
calibration to 5%. We note that both the astrometric and photometric precision increase as 
the square root of the ratio of the satellite PSF intensity to the speckle noise background, 
i.e. Fa/^raask if both the PSF and mask are dominated by the second order halo term at 
the satellite PSF separation. Since this is the dependance that was found in sect. 4 for the 
unaberrated case, it is tempting to conclude that the second term of Eq. 11 introduces both 
astrometric and photometric errors, as was expected from the fact the Fa is complex. 

A satellite PSF relative intensity of ~ 100 times the background speckle noise seems 
to be a good choice to avoid over-luminous satellite PSFs while offering good registering 
and photometric accuracy (0.3 pixel RMS (~1/15 \/D) and 5%). Such an astrometric 
precision (0.003 arcsec for diffraction limited H-band images acquired with a 10-m telescope) 
is sufficient to conffim proper motion of nearby systems (typically 0.01 to 0.1 arcsec/year) 
in one year.^ Due to detector dynamic range limits, to minimize satellite PSF intensities, 
it might be more convenient to use a large number of fainter satellite PSFs to average 
the speckle noise bias, albeit at the cost of a greater FOV contamination, or choose a mask 
frequency to move the satellite PSFs to wide separations, where the speckle noise is generally 
smaller. 



^For speckle noise limited companion detections, the derived astrometric precision from Sect. 5 simulations 
also mean that companions need to be at least 100 times brighter than the background speckle noise to 
obtain 1/15 \/ D center determination accuracy and allow their proper motion follow-up in one year for 
typical nearby systems - fainter companions will require a bigger time interval or longer integration time, 
assuming that the object S/N is increasing with integration time. 
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6. Laboratory Experiment 

A transmissive amplitude grating lias been tested on the Universite de Montreal high 
contrast imaging testbed. Commercial lenses were used to make a f/31 imaging system, 
yielding, at 1.625 /xm and with a 3% filter bandpass, a PSF having 2.7 pixels per \/ D. The 
pupil diameter was 13.1 mm. A transmissive amplitude grating mask having 150 ± 10 /im 
thick wires and 595 ± 15 //m spacing was introduced 10 mm from the pupil plane. This slight 
offset was necessary due to the mask and pupil stop mounts. The geometric pupil diameter 
at the grating mask conjugated plane is calculated to be 12.75 mm. With this setup, the 
diffraction envelope of a single wire drops to 0.97 of its peak at a separation of 17.1 A/-D, 
the location of the first satellite PSF pair. Neglecting the satellite PSFs chromaticity, it is 
less than 0.6 A/D, it is thus expected that this grating will produce satellite PSFs having a 
relative intensity of 0.97 x (150/595)^ = 0.062 ± 0.012 at 17.1 ± 0.5 \/D separation. The 
pinhole was moved and an aberration phase plate was introduced in the beam, producing a 
Strehl ratio 25% lower, to confirm that these satellite PSFs do track the PSF core position 
and intensity. All acquired images were unsaturated and reduced by subtracting a dark 
frame and dividing by a fiat field. The two brightest satellite PSFs are 160 times brighter 
than the background noise. Fig. 5 shows the resulting PSFs. Satellite PSFs have a relative 
intensity of 0.058 and are situated at 16.9 \/ D from the PSF center, both values being 
consistent to the calculated values to la accuracy, and they track the PSF core position and 
its intensity, in agreement with what was expected. In all above experiments, it was possible 
to register the PSF to less than 0.1 pixel precision (1/30 \/D) with the average position 
of two symmetric satellite PSFs and the relative intensities are constant to several percent 
(see Table 4), all consistent with the expected accuracy for satellite PSFs that are 100-1000 
times the background noise (see Table 3). 



7. On-Sky Implementation 

If we consider the Gemini Altair AO system (Herriot et al. 1998) and the NIRI near- 
infrared camera (Hodapp et al. 2000) for a possible implementation, satellite PSFs ~ 10 
magnitudes fainter than the central peak should be ~100 times brighter than the speckle 
noise at 4 arcsec (Marois et al. 2006). If we put a transmission amplitude grating just in 
front of Altair's deformable mirror, where the wavefront diameter is 84 mm, a wire spacing 
and thickness of 2 mm and 18 /im would respectively be required to produce satellite PSFs 
10 magnitudes fainter at 1.7 arcsec intervals. Such a configuration has 99.1% transmission 
(see Fig. 6 for a simulated example). 

For high contrast imaging seeking for detection limits of the order of 20 magnitudes 
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at less than 1 arcsec, e.g the Gemini Planet Imager (GPI) and the VLT Planet Finder 
(VLTPF) projects, wire thickness will become an issue. A potential solution is to increase 
the wire spacing to place the satellite PSFs at 5 \/ D. Using the same 18 /xm thick wires 
and wavefront diameter as the above example, it would produce satellite PSFs that are 15 
magnitudes fainter than the primary. For more advanced coronagraphic projects aiming for 
Earth- like exoplanet detections at 24 magnitudes contrast, a simple approach is to put the 
wires directly on the primary mirror. A series of ten 18 microns thick wires across an 8-m 
diameter telescope will produce satellite PSFs 23.2 magnitudes fainter than the primary at 
10 \/D (99.998% transmission). 

Another solution is to generate satellite PSFs using a small periodic phase aberration. 
In a conventional adaptive optics system with a separate wavefront sensor this could be 
achieved by adding a sinusoidal offset to the control point that the wavefront sensor at- 
tempts to achieve. Most AO systems operate with their wavefront control loop seeking not a 
flat wavefront in the sensor but a wavefront offset by the conjugate of any non-common-path 
aberrations between the wavefront sensor and science channel. Adding an additional sinu- 
soidal component to this offset would result in a time-averaged PSF containing the satellite 
PSFs. As long as the sinusoid is lower in amplitude than the typical residual wavefront, 200- 
300 nm for current AO systems (van Dam & Macintosh 2003) and 50-100 nm for extreme 
AO (Macintosh et al. 2004; Mouillet et al. 2004), and the wavefront sensor is operating in a 
linear regime this should have little effect on other components of the image. For systems 
without a separate wavefront sensor - e.g. the Terrestrial Planet Finder Coronagraph, which 
intends to use science-camera images on a reference target to control the deformable mirror 
via iterative phase diversity, the satellite images could be built into the target wavefront or 
focal-plane image that the phase diversity image sequence converges to. The satellite PSFs 
can then be used to estimate the PSF changes between the reference and science targets. 



8. Discussion 

Simulations presented in this paper are monochromatic. For bandpasses of finite size, 
off-axis satellite PSFs will be elongated and thus fainter relative to the core. The satellite 
PSF elongation is simply proportional to the filter bandpass AA/A times the satellite PSF 
angular separation. Additionally, for ground-based imaging, if no atmospheric dispersion 
corrector is used, the elongation of the satellite PSFs will not necessarily points towards 
the PSF center due to differential atmospheric refraction between the blue and red side of 
the bandpass (Filippenko A. V. 1982; Roe H. G. 2002). For H-band observations within 
±2HA, the amplitude of this effect can be several tens of mas, i.e. of the order of one 
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A/-D, as observed by precise registration of two near-infrared narrow-band PSFs made with 
the TRIDENT camera at CFHT (Marois et al. 2005). However, satelhte PSF chromaticity 
and differential atmospheric refraction do not influence the astrometry precision since the 
midpoint of the line connecting the centers of opposite satellite PSFs still coincides with 
the PSF center. For large bandpasses and/or satellite PSFs at large angular separations, 
the satellite PSF photometry can be partially correlated with a Strehl ratio variation since 
the required aperture necessary to estimate the satellite PSF intensity will need to be larger 
than a \/ D diameter; it would thus be better to choose a satellite PSF at small separation 
to minimize this effect. 

Distortions coming from optics after the added mask can change the relative position 
of the satellite PSFs. The satellite PSF position will thus need to be calibrated using an 
internal source or on-sky observations. To obtain precise relative astrometry of detected 
companions, distortion corrections for the entire optical system is also required. For ground- 
based imaging, two effects are potentially not correctable with the satellite PSF technique. If 
a wide bandpass is used, spectral differences between the companion and the primary coupled 
with differential atmospheric refraction can introduce small astrometric errors as a function of 
airmass. This effect is negligible for IFUs or narrowband imaging and can, for cases when it is 
important, be accounted for by analysis the companion spectrum. Anisoplanatism variations 
can also degrade the companion Strehl ratio without affecting that of the primary, thus 
possibly introducing companion photometric errors. Since high contrast imaging projects 
usually have small FOV (several arcsec wide) compared to standard AO anisoplanatism 
angle (~30 arcsec radius), such effect is relatively small. 

The satellite PSFs will occupy a small portion of the FOV. Once the PSF is registered 
and flux normalized, they could be removed to avoid masking faint sources. Since the 
angular separation and size of the satellite PSFs is chromatic and fixed with respect to the 
instrument, a simple solution to remove them while conserving the flux of nearby sources 
would be to use the simultaneous spectral differential imaging SSDI (Racine et al. 1999; 
Marois et al. 2000; Sparks & Ford 2002; Biller et al. 2004; Marois et al. 2005) and/or the 
angular differential imaging techniques ADI (Marois 2004; Liu 2004; Marois et al. 2006). If 
such techniques still leave important localized residuals, the satellite PSFs can simply be 
masked. 



9. Conclusion 

It was shown both analytically and numerically that periodic amplitude and phase masks 
conjugated to the pupil plane with and without a coronagraph can be used to produce 
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satellite PSFs that track the PSF core position and its intensity. The astrometric and 
photometric precision is proportional to the square root of the ratio of the satelhte PSF 
intensity to the background speckle noise. Satellite PSFs 100 times brighter than the speckle 
noise achieve a registering and photometric accuracy of 0.3 pixel RMS (~1/15 \/ D) and 5% 
respectively. Such a precision is sufficient to properly combine images and to confirm proper 
motion of nearby systems within one year with AO on a 10-m telescope. 

As illustrated by the laboratory experiment, this technique can be easily implemented 
by using a transmissive amplitude grating. The experiment has validated that an amplitude 
mask conjugated to the pupil plane produces such satellite PSFs and that they track the 
PSF core position and intensity while demonstrating the simplicity of the method. 

This technique could play an important role in current AO systems and next generation 
extreme adaptive optics coronagraph systems like the Gemini Planet Imager (Macintosh et 
al. 2004, 2006) or the VLT Planet Finder (Mouillet et al. 2004) as well as in future high 
contrast coronagraphic space observatories. 
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Sciences and Engineering Research Council, Canada and from the Fonds Quebecois de la 
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Table 1: Astro metric and Photometric Precision with Mask Offset {\/ D = 4 pixels, no 
aberration). 



Maslc type 


Offset (period) 


Coro. 


Axe 


Aye 


-^satel/-^PSF 


Phase 





no 


-0.31 


0.00 


0.00100 




1/4 


no 


0.00 


0.00 


0.00099 







yes 


0.00 


0.00 


0.00100 




1/4 


yes 


0.00 


0.00 


0.00099 


Amphtude 





no 


0.00 


0.00 


0.00100 




1/4 


no 


0.00 


0.00 


0.00088 







yes 


0.00 


0.00 


0.00100 




1/4 


yes 


0.00 


0.00 


0.00100 


Grating 





no 


0.00 


0.00 


0.01211 




1/4 


no 


0.00 


0.00 


0.01237 







yes 


0.00 


0.00 


0.01240 




1/4 


yes 


0.00 


0.00 


0.01240 
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Table 2: Simulated Astrometric and Photometric Precision {\/ D = 4 pixels). 



Aberration 


Coro. 


With 


AXC rms 


AYC rms 


-^satelZ-^PSF 


-^satclZ-^PSF 


-^satclZ-^PSF 






150nm RMS 


(pixel) 


(pixel) 




rms (%) 


expected 


Phase 


no 


no 


- 


- 


0.00100 


- 


0.00100 


Phase 


no 


yes 


0.19 


0.18 


0.00100 


0.00001 (1.0) 


0.00100 


Phase 


yes 


yes 


0.15 


0.10 


0.00099 


0.00001 (1.0) 


0.00100 


Amplitude 


no 


no 


- 


- 


0.00100 


- 


0.00100 


Amplitude 


no 


yes 


0.24 


0.20 


0.00102 


0.00005 (5) 


0.00100 


Amphtude 


yes 


yes 


0.19 


0.18 


0.00103 


0.00005 (5) 


0.00100 


Grid 


no 


no 


- 


- 


0.01240 


- 


0.01235 


Grid 


no 


yes 


0.12 


0.15 


0.01198 


0.00010 (0.8) 


0.01235 


Grid 


yes 


yes 


0.05 


0.03 


0.01235 


0.00010 (0.8) 


0.01235 
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Table 3: Satellite PSF Intensity vs Astrometric and Photometric Precision {\/D = 4 pixels). 



Satellite PSF int. 
to speckle noise ratio'' 



Coro. AXC rms AYC rms /satci/-^PSF 



-^satclZ-^PSF 

rms (%) 



10 

16 

100 

160 

1,000 

1,600 

10,000 

16,000 



no 


0.78 


yes 


0.75 


no 


0.24 


yes 


0.19 


no 


0.07 


yes 


0.06 


no 


0.02 


yes 


0.02 



0.93 


0.00013 


0.00002 (15) 


0.80 


0.00011 


0.00002 (18) 


0.20 


0.00102 


0.00005 (5) 


0.18 


0.00103 


0.00005 (5) 


0.06 


0.01003 


0.00018 (1.8) 


0.05 


0.01005 


0.00014 (1.4) 


0.02 


0.10214 


0.00092 (0.9) 


0.02 


0.10156 


0.00067 (0.7) 



^ In all cases, satellite PSF intensities are 0.1% fainter than the PSF peak intensity. The satellite PSF 
intensity to speckle noise ratios are higher when using a coronagraph since, for these cases, the speckle noise 
is attenuated by the coronagraph by a factor 1.6 in our simulations. 



Table 4: Laboratory Experiment. 
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PSF 


Satellite PSFs 










XC YC 


XC YC 


AXC 


AYC 


< -^satel/-^PSF > 


Unaberrated 


199.93 199.98 


199.90 200.05 


0.03 


-0.07 


0.055 


Aberrated 


200.08 199.69 


199.99 199.68 


0.09 


0.01 


0.061 


Shifted 


200.45 209.04 


200.45 209.07 


0.00 


-0.03 


0.059 
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Fig. 1. — Central plot: Simulated monochromatic PSF with added sine wave phase mask 
(solid line). The dashed line shows the same simulation with an added tilt to show that the 
satellite PSFs track the PSF core position. Top panel and bottom plot: the same simulation 
with an added 150 nm RMS phase aberration. Dotted lines show the satellite PSFs expected 
positions and intensities. Top panel is shown on a logarithmic scale. 
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Fig. 2. — Same as Fig. 1 for an added amplitude sine wave mask. 
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Fig. 3. — Same as Fig. 1 for an added transmissive amplitude grating. This figure shows the 
first two bright satelhte PSFs. Numerous fainter satelhte PSFs exist at wider separation. 



-23- 



'm 

c 



N 



1^ 

o 



-t-J 

en 



N 




-30 



-20 



-10 10 

Ang.ular Separation (A/D) 



.20 



30 



Fig. 4. — Same as Fig. 1 for a Gaussian apodized pupil coronagraph with a transmissive 
amplitude grating. 
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Fig. 5. — Same as Fig. 1 for the laboratory experiment using a transmissive amplitude 
grating. 
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Fig. 6. — Observed 30 s CH4-band Gemini Altair/NIRI saturated PSF (left, star has mn = 
5) and with simulated satellite PSFs (right). Simulated satellite PSFs are for a series of 42 
wires having 18 fira thickness across a 84 mm pupil. Satellite PSFs are ~ 10 magnitudes 
fainter than the primary and are located every ~1.7 arcsec interval along the X axis. They 
are produced using non-saturated data acquired in the same filter and normalized to the 
same integration time as the saturated image. Satellite PSFs have approximately the same 
intensity contrast in the FOV since the sine modulation FWHM from the wire thickness (see 
Eq. 13) is ~ 9 times the FOV (FOV is 22.4 x 22.4 arcsec). The satellite PSF elongation is 
simulated for the filter 6.5% bandpass. PSF is saturated inside 0.8 arcsec radius. 



